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Robert  E.  L.  Turner 

Internal  Waves  in  Stratified  Atmospheres 

Work  was  continued  on  the  study  of  internal  waves  in  stratified  fluids  using  both  an¬ 
alytical  and  computational  techniques.  Attention  was  particularly  directed  to  waves  in  a 
two-fluid  system,  confined  between  horizontal  walls  and  acted  on  by  gravity.  The  underly¬ 
ing  equations  are  the  Euler  equations  for  nondiffusive  flov/.  For  each  base  configuration  of 
flmd  depths  and  densities  it  was  shown  by  Amick  and  Turner  in  [1]  that  there  axe  global 
branches  of  solitary  wave  solutions  starting  with  infinitesimal  waves  corresponding  to  a 
critical  speed  for  the  system.  Open  questions  regarding  the  limiting  form  of  these  waves 
led  to  various  projects  [2], [3], [4]  which  were  continued  or  completed  during  the  grant  period. 

Amick  and  Turner  conducted  analytical  investigation  of  all  small  waves  in  the  two  fluid 
system,  particularly  surges.  A  surge  is  a  flow  which  approaches  different  limiting  parallel 
flows  as  the  horizontal  coordinate  approaches  infinity  to  the  left  and  to  the  right.  No 
rigorous  analytical  studies  of  surges  had  been  presented  in  the  literature  before,  and  now  in 

[2]  we  have  a  complete  description  of  all  small  amplitude  interfaces  between  two  irrotational 
regimes  of  differing  densities  under  the  influence  of  gravity.  We  used  a  dyneunical  systems 
approach  and  a  center  manifold  reduction  to  model  the  behavior  of  small  waves.  In  this  way 
we  could  parametrise  all  ’small’  solutions  of  the  full  system,  a  quasi-linear  elliptic  problem 
arising  from  the  Euler  equations  through  suitable  changes  of  coordinates,  by  means  of 
solutions  of  an  ordinary  differential  equation  of  second  order  in  a  center  manifold.  In  fact, 
the  existence  of  an  invariant  for  the  flow  allows  a  final  reduction  to  a  first  order  differential 
equation,  the  heteroclinic  orbits  for  which  correspond  to  surges.  The  work  for  this  project 
led  to  the  broader  study  of  functional  equations  of  the  type  arising  in  hydrodynamics.  In 

[3]  we  investigate  mappings  in  weighted  Holder  spaces  which  are  useful  for  the  analysis  of  a 
general  class  of  functioned  equations.  These  equations  include  those  arising  in  the  internal 
wave  problem  as  well  as  those  from  surface  wave  problems  with  surface  tension  [5]  and,  we 
expect,  will  be  useful  in  other  wave  problems. 

A  project  completed  under  the  grant  was  a  numerical  study  [4],  done  jointly  with 
Jean-Marc  Vanden-Broeck,  to  determine  the  limiting  forms  possible  for  the  waves  shown 
to  exist  in  the  paper  [1].  From  that  paper  we  knew  that  along  a  branch  of  solutions  the 
solitary  waves  would  have  to  broaden  indefinitely  or  the  slopes  of  the  streamlines  in  the 
flow  would  have  to  steepen  approaching  the  vertical.  The  problem  was  recast  in  terms  of 
linked,  nonlinear  integrodifferential  equations  in  one  variable  parametrising  the  interface. 
The  equatioiiS  were  discretized  and  the  resvdting  nonlinear  equations  solved  using  the 
Newton-Raphson  method.  It  was  shown  that  for  certedn  ranges  of  the  density  and  depth 
parameters  the  solitary  waves  appeared  to  broaden  indefinitely.  For  example,  with  density 
ratios  like  that  of  mr  and  water  and  with  equal  fluid  depths,  broad  waves  were  seen.  This 
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was  evidence  for  the  existence  of  stirges,  a  surge  being  roughly  ’half’  of  an  extremely  broad 
wave. 
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Summary  of  Scientific  Progress 

Robert  C.  Rogers.  Postdoctoral  Associate 


1  Self-Effect  Problems  in  Electro-Magneto-Elasticity 

Professor  Rogers’  v  irk  in  [6]  and  with  S.S.  Antman  in  [l]  is  part  of  a  program  to  extend  the  modern 
existence  and  regularity  theory  of  elasticity  to  problems  in  steady-state  electromagnetism.  These  works 
examine  problems  involving  polarizable,  magnetizable,  conducting,  deformable  materials  with  monotone 
electromagnetic  response.  The  constitutive  equations  are  assumed  to  have  the  general  form: 

T(x)  =  T(F(x),e(x),h(x).x), 

p(x)  =  p(F(x),  e(x).  h(x),  x),  (1.1) 

in(x)  =  rii(F(x),e(x),h(x),x), 
jW  =  j(F(x),  e{x),  h(x),  x). 


Here  T,  p,  m,  and  j  are  the  Second  Piola-Kirchhoff  stress  and  the  Lagrangian  forms  of  the  polarization, 
magnetization,  and  electric  current  respectively;  and  F,  e,  and  h  are  the  deformation  gradient  and  the 
Lagrangian  versions  of  the  electric  and  magnetic  fields  respectively.  These  coupled,  nonlinear  relations  allow 
for  such  effects  as  electrostriction  and  magnetostriction. 

A  suitable  form  of  the  Biot-Savart  law  implies  that  the  magnetic  field  generated  by  the  conducting, 
magnetized  body  satisfies  the  following  integral  equation: 


h(*)  = 


j(F(x),e(x),h(x) ,  x)(g  -  x) 


In  [l]  it  was  shown  that  under  appropriate  hypothescson  j,  h  has  the  representation 


h(x)  =  Vt6(x)-f«(F,e,VV>,x), 


(1.3) 


where  V*  is  called  the  magnetic  potential  and  (P,e,  V^)  i— »  K(F,e,  is  compact  on  appropriate 

spaces.  Composing  (1.3)  with  the  constitutive  equations  (1.1),  one  can  write  the  balance  laws  of  elasticity 
and  electro-magneto-statics  in  the  form: 

DivT(F(),V0(-),VV.(-),F(x),V^(x),VV>(x),x)  =  -f(x), 

Div  d(F(  ),  V^(  ),Vv;(  ),F(x),\'^(x),VVi(x),x)  =  <r(x),  (1.4) 

Divb(F(),V^(),VV-(),F(x),r<J(x).Vit(x),x)  =  0. 


Here  (t>  is  the  electric  potential,  and  the  definitions  of  d  and  b  and  f  can  be  formed  from  e,  h,  and  F  in  the 
usual  way. 

The  existence  of  a  large  class  of  solutions  for  the  system  (1.4)  was  shown  in  [l]  using  techniques  of 
pseudo-monotone  operators  (c/.  [ll])  The  key  observation  is  that  the  constitutive  equations  T,  d,  and  b  are 
compact  in  the  global  variables  and  monotone  in  the  local  variables.  In  [6]  a  more  general  class  of  boundary 
conditions  was  considered,  but  greater  restrictions  were  put  on  the  constitutive  equations. 


2  Modeling  of  Ferromagnetic  Materials 

Professor  Rogers  has  made  several  advances  in  the  theory  of  ferromagnetic  materials.  In  [6]  the  classi¬ 
cal  micromagnetic  model  was  investigated,  and  an  existence  theorem  was  proved  for  coupled,  unshielded, 
magneto-elastic  materials  under  very  general  boundary  conditions.  The  proof  of  the  theorem  involves  meth¬ 
ods  of  convex  analysis  and  weak  ci  niinuity.  Nonlocal,  shape-dependent  effects  are  taken  into  account. 

In  addition,  new  models  of  ferromagnetic  materials  have  been  formulated  and  studied.  In  ^2i  a  nonlocal 
constitutive  relation  is  considered,  and  an  existence  theorem  is  proved  using  pseudo-monotone  operator 
methods. 
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In  .8]  a  nonlocal  energy  model  is  studied.  The  new  model  allows  for  the  consideration  of  measure-valued 
magnetizations,  which  are  suggested  by  recent  studies  of  Erickses  [15-,  Chipot  L  Kinderlehrer  [12], 
and  Ball  4c  James  [9]  on  twinning  in  crystals.  The  energy  of  a  ferromagnetic  body  with  magnetization  m 
in  an  applied  field  ho  is  given  by 


{W{m)  -  m -ho}*!!' -i- x(>^)-  (2-1) 

The  first  integral  is  the  field  energy,  h,(m)  is  the  resultant  magnetic  field  induced  by  the  magnetization  m. 
In  the  second  integral,  the  first  term  is  the  magnetization  energy  density.  We  assume  that  W  is  nonconvex. 
In  an  isotropic  material  we  expect  W  to  be  minimized  for  some  prefered  nonzero  modulus  of  magnetization 
A<o-  The  second  term  in  the  integral  is  the  interaction  energy, and  the  final  term  x  is  the  exchange  energy. 
In  the  theory  of  micromagnetics  an  exchange  energy  of  the  form 


Xt(m)  =  c  /  lVm(x)l*du, 

is  used.  In  [8]  the  exchange  energy  is  modeled  by  a  nonlocal  term  of  the  form 


Xj(m)  =  C(7) 


Ul/, 


(m(x)  - 


|x-y| 


dv^ 


(2.2) 


(2.3) 


The  new  model  for  the  exchange  energy  is  designed  to  penalize  local  deviations  in  the  magnetization 
while  allowing  for  consideration  of  a  wider  class  of  magnetizations  than  that  implied  by  an  energy  of  the 
form  (2.2).  In  particular,  one  can  consider  measure- valued  distributions  of  magnetization.  The  use  of  such 
distributions  is  motivated  by  the  work  on  twinning  in  crystals  cited  above. 

If  no  exchange  energy  is  present,  the  following  oscillating  sequence  of  magnetizations  minimizes  the  energy 
absolutely  for  -Mo/3  <  'Ho  <  Mo/3: 


X  lv 


=  Mo/,1 ,32jfl..(—:-)k,  *  £  J?i(0), 

V  a  ^  avMo  *  J 


(2.4) 


Here  i  is  a  unit  vector  perpendicular  to  k  and  /«  is  a  function  of  period  1  defined  on  the  real  line  for  0 
by 


0  <  *  < 

0  <  z  <  1. 


<  ® 
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The  sequence  does  not  converge  strongly  in  i*  but  does  converge  weakly  to  a  uniform  “average”  mag¬ 
netization: 


Note  that  the  average  magnetization  does  not  solve  the  original  problem.  (It  solves  what  is  known  as  the 
relaxed  variational  problem.)  The  idea  introduced  in  the  references  on  twinning  cited  above,  is  that  one  can 
think  of  the  Young’s  measure  of  the  minimizing  sequence  as  solving  the  problem.  Essentially,  the  Young’s 
measure  describes  the  oscillations  of  a  weakly  converging  sequence.  For  instance,  the  Young’s  measure 
associated  with  the  sequence  m-'  would  be 


t/ 


-Mok). 


Here  5(v)  is  a  Dirac  measure  centered  at  v.  The  measure  i/  reflects  the  probability  finding  the  values  of  the 
sequence  in  either  of  the  states  -Mok  in  any  small  ball  in  Bi(0).  The  center  of  mass  of  the  measure  is 
equal  to  the  weak  limit  of  the  sequence,  the  average  magnetization. 

To  see  how  the  exchange  energy  xz  penalizes  oscillations  we  examine  how  it  is  affected  by  the  sequence 
m*.  Some  routine  calculations  give  us  the  result 


lim  xj(ni^)  =  A' 

J  — •  00 


|^)(Mo^Ho/3)*j, 


(2.5) 
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where  K  is  a  constant  dependent  on  7.  Thus,  xi  penalircs  the  variance  of  the  Young’s  measure. 

The  point  of  view  taken  in  the  studies  of  twinning  mentioned  above  that  Young’s  measures  can  be 
thought  of  as  solutions  of  the  minimization  problem.  (This  point  of  view  descends  directly  from  the  ideas  of 
generalized  curves  of  L.C.  Young.)  In  [8j  it  is  shown  that  the  introduction  of  the  nonlocal  exchange  energy- 
leads  to  multiple  measure-valued  minima  of  the  energy. 

In  [7]  the  following  model  problem  was  examined  order  to  get  a  further  idea  of  the  effects  of  penalizing 
the  variance  of  a  measure- valued  solution.  Let  u  =  (ui, . . . ,  «n)  be  and  element  of  R".  The  following  energy 
was  minimized 


Here  W  is  the  nonconvex  energy 


yv(u) 


I  ^(n-uo)*  t»>0, 

\  ^{v  +  uoY  i)<0, 

y(u-  sgn(t))uo)’, 


(2.6) 


(2.7) 

(2.8) 


and  uo,7<o,  Qi,a2,ot3,<»4  are  positive  constants.  The  first  term  is  analogous  to  the  magnetization  energy,  the 
second  to  the  field  energy,  the  third  to  nonlocal  exchange  energy,  and  the  fourth  to  the  interaction  energy 
with  Tio  representing  the  strength  of  the  applied  field.  Figure  1  represents  the  minimizers  for  n  =  8.  We 


Figure  1:  Solution  curves  for  u  G  R*.  Solution  segments  lie  between  the  dashed  lines 
p  =  (o^Tfo  ±  oi«o)/(Q2  -  “s)-  All  solution  segments  are  parallel.  The  upper  and  lower  half  lines  of 
solutions  are  on  the  lines  fi  =  {a^Tio  x:  aiUo)/(oi  -r-  oj). 

have  graphed  u  =  ’*•  ’'‘’hich  we  think  of  as  analogous  to  the  center  of  mass  of  the  Y'oung’s  measure  as  a 

function  of  Rq. 

Raising  Q3  (the  coefficient  of  the  moment  term,  out  model  for  the  exchange  energy)  causes  the  boundary- 
lines  to  become  more  vertical.  If  we  take  the  point  of  view  that  as  Ro  is  slowly  changed  we  will  stay  on 
a  solution  curve  until  forced  to  change  by  contact  with  the  boundary,  a  more  vertiezd  boundary-  causes  the 
inner  loops  of  solutions  to  be  much  less  stable  than  the  outermost  loop,  thus  making  the  solutions  more  like 
the  hysterysis  diagram  of  a  hard  ferromagnetic  material. 
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3  Problems  in  Viscoelasticity 

Recent  work  of  Professor  Rogers  with  J.A.  NOHEL  and  A.E.  Tzavaras  [5’  initiated  a  program  of  investi¬ 
gating  the  existence  of  global  weak  solutions,  in  the  class  of  bounded  measurable  functions,  for  the  Cauchy 
problem  for 

“  "*’  1 6  R,  t  >  0,  (3J) 

1),  =  o-,. 


I  6  R,  t  >  0, 


with  initial  conditions 


t£;(*,0)  =  too(*),  t;(x,  0)  =  t)o(a:),  xER.  (3.2) 

Here  the  function  <r{x,t)  is  determined  by  the  history  of  through  the  constitutive  assumption 


'■(*■0  =  +  [  Hi- T)rl){w{x,T))dt 

Jo 


The  given  functions  v’(u;),V’(w)  and  k(t)  are  assumed  to  be  smooth  and,  in  addition, 

Vi'(ttj)  >0,  to  €  R,  (3.4) 

so  that  the  structure  of  (3.1)  is  hyperbolic.  This  equation  is  a  widely  accepted  model  for  important  phe¬ 
nomena  in  viscoelastic  solids. 

In  [5]  the  important  special  case  il>  =  <p,  is  considered  when  the  data  lOo,  wo  6  L°°(R)  n  i*(R).  The 
problem  becomes 

~  \  ^  r*U4  \  i  (  \\  A  *  G  R,  1  >  0,  (3.5) 

t»,  =  v>(ti)), -l-/oA:(t-T)v>(io(-,T)),dr, 

io(*,0)  =  too(x),  t;(i,0)  =  uo(»).  X  G  R. 

The  constitutive  function  tp  is  assumed  to  satisfy 

^  :  R  — >  R  is  a  twice  continuously  differentiable 
^  function  such  that  v’*(«')  >  0i  E  R;  gx 

if  has  a  single  inflection  point  &t  w  =  Wi  and  is  \  •  f 

convex  on  (roi,oo)  and  concave  on  (— oo,u;j). 


The  kernel  k  satisfies 


and  the  data  ioo(x),  uo(x)  satisfy 


jfe;[0,oo)-R,  jfeeC*[0,oo), 
i«o(3:),  vo{x)E  L^{R)nL\R). 


The  following  theorem  was  proved. 

Theorem  3.1  Let  the  hypotheses  (3.6)-(3.8)  be  satisfied.  Given  T  >  0,  there  exists  a  weak  solution 
{«;(x,  t),  i’(x,  l)}  of  (3.5)  on  R  X  [O.Tj,  such  that 

(w,v)  E  I*([0,T];X*(R))nI“(R  x  [0,T]).  (3.9) 

DiPerna  [13'j  established  similar  results  in  the  case  of  1-d  elasticity  and  1-d  isentropic  gas  dynamics. 


4  Problems  in  Compensated  Compactness  and  Weak  Continuity 

The  techniques  of  compensated  compactness  and  the  characterization  of  weakly  continuous  functionals  under 
various  differential  constraints  has  been  an  important  tool  in  some  recent  developments  in  partial  differential 
equations.  In  particular,  the  Div-Curl  Lemma  was  instrumental  in  the  work  of  Tartar  [16  and  DiPerna 
[13]  on  conservation  laws  and  in  our  work  on  weak  solutions  in  viscoelasticity  [5_;  and  the  characterization 
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of  the  null  Lagrangians  in  the  setting  of  the  calculus  of  variations  (cf.  Edelen  [14  )  was  central  to  the  work 
of  Ball  [10]  on  existence  theorems  in  three-dimensional  elasto-slatic  polyconvex  functions. 

Professor  Rogers,  in  work  with  J.  Robbin  and  B.  Temple  |3'  devised  a  method  sufficient  for  construct¬ 
ing  weakly  continuous  polynomials  via  wedge  products.  This  method  produces  polynomials  of  degree  as 
high  as  the  domain  space,  the  maximum  degree  possible.  With  this  method  one  can  reproduce  all  of  the 
weakly  continuous  functions  in  the  div-curl  case,  the  case  of  Maxwell's  equations  of  electrodynamics,  and  the 
variational  case.  As  Tartar  indicates  jl7],  the  Quadratic  Theorem  can  be  used  to  prove  such  a  theorem  when 
£*  techniques  apply.  The  proof  in  [3]  is  based  on  the  idea  of  decomposing  weakly  convergent  sequences  into 
weakly  convergent  and  strongly  convergent  parts  using  a  version  of  the  Hodge  decomposition.  The  result  is 
more  general  than  that  obtained  directly  from  the  £*  theory.  In  addition,  the  proof  gives  additional  physical 
insight  into  the  phenomenon  of  weak  continuity  and  the  interaction  of  oscillations. 

The  work  on  wedge  products  led  to  a  more  general  sufficient  condition  for  the  weak  continuity  of  poly¬ 
nomials.  The  condition  (given  in  [4])  relates  weak  continuity  to  the  wave  cone  and  characteristic  set  of 
the  differential  constraints.  It  is  both  necessary  and  sufficient  when  the  polynomial  is  quadratic  or  when 
a  certain  orthogonality  condition  is  satisfied.  Since  the  condition  is  more  general  than  the  wedge  product 
condition  it  applys  to  all  of  the  classical  examples  of  weakly  continuous  polynomials. 

This  condition  clarifies  the  relationship  between  the  directions  in  the  wave  cone  (associated  with  the 
dependent  variables  of  the  weakly  convergent  sequence)  and  the  characteristic  directions  (associated  with 
the  independent  variables  of  the  sequence).  The  Legendre- Hadamard  condition  involves  only  the  wave  cone, 
but  for  degree  greater  than  two,  the  coupling  of  these  directions  to  the  characteristic  directions  is  crucial. 

The  method  of  proof  for  the  general  condition  is  to  obtain  conditions  under  which  one  can  iterate  the 
quadratic  theorem  of  Tartar.  In  order  to  do  this  we  characterize  fully  the  first-order  differential  constrmnts 
induced  in  quadratics  by  first-order  differential  constraint  on  their  factors. 
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Peristaltic  Pumping  with  Heat  Release 
and  Heat  Conduction 

M.  C.  Shen  and  Oalin  Tang 


Peristalsis  was  originally  a  physiological  pumping  mechanism  which 
transports  liquid  along  a  duct  by  a  large-amplitude,  wavelike  motion  of  the 
duct  walls;  the  esophagus,  urethelia  and  urethra  work  that  way.  It  was  picked 
up  by  medical  engienering  for  antiseptic  liquid  transport  and  by  nuclear 
engineering,  because  it  allows  radioactive  liquid  to  be  transported  without 
corrosion  of  the  pumping  machinery.  In  the  latter  context,  in  particular, 
radioactive  heat  release  and  heat  transfer  become  important  and  little 
knowledge  on  their  influence  had  oeen  available  before  Dr.  Tang  took  up  the 
subject  in  his  Ph.D.  research  under  the  supervision  of  Professor  M.  C.  Shen. 

For  the  elucidation  and  prediction  of  such  peristalsis.  Profess  Shen  and 
Dr.  Tang  have  developed  two  eymptotic  methods  which  are  fully  nonlinear  as 
needed  in  applications;  one  is  designed  for  small  Reynolds  numbers  and  the 
other,  for  waves  long  compared  to  the  duct  diameter  over  the  whole  range  of 
Reynolds  number  of  practical  interest. 

In  order  to  demonstrate  the  validity  of  these  approximations  conclusively, 
they  have  also  established  the  rigorous  theory  of  the  nonlinear  Overbeck- 
Boussinesq  equations  governing  peristaltic  fluid  motion  with  heat  conduction 
and  heat  release  on  a  wide  basis.  By  the  help  of  the  large  apparatus  of 
inequalities  needed  to  make  either  degree  theory  or  Galericin  limits  work,  they 
have  proved  the  wel Iposedness  of  not  only  the  basic  periodic  pumping  state, 
but  also  of  the  problems  for  unsteady  and  nonperiodic  motions  of  the  duct 
walls.  These  theorems,  moreover,  have  permitted  them  to  prove  the  stability  of 
the  basic  pumping  state  by  energy  estimates.  It  appears  hard  to  avoid  the 
conclusion  that  this  research  will  stand  as  the  definitive  theoretical  work  on 
peristalsis  with  heat  release  and  heat  conduction. 

In  addition.  Dr.  Tang  has  employed  the  mathematical  knowledge  gained  to 
compute  both  exact  and  approximate  solutions  of  the  Overbeck-Boussinesq 
equations  and  has  shown  thereby  that  his  approximations  have  satisfactory 
accuracy  over  the  whole  range  of  practical  interest.  Any  nuclear  or  medical 
engineering  project  can  therefore  use  his  asymptotic  methods  confidently  to 
predict  the  parameter  ranges  over  which  peristaltic  pumping  will  work  without 
back-flow  or  significant  pockets  of  trapped  liquid.  A  number  of  mathematical 
papers  for  the  publication  of  all  this  work  are  in  preparation. 
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Research  by  S.  M.  Sun  under  the  supervision  of  M.  C.  Shen 


1.  Viscous  magnetic  fluid  flow  down  an  inclined  plane. 

The  problem  deals  with  a  viscous  magnetic  fluid  flow  down  an  inclined 
plane  under  the  influence  of  an  applied  magnetic  field  parallel  to  the  plane. 
The  commonly  used  method  is  to  find  a  critical  Reynolds  number  for  stability 
in  the  long  wave  limit  via  linearized  equations.  Our  approach,  instead,  leads 
to  a  nonlinear  evolution  equation  governing  the  development  of  a  weakly  non¬ 
linear  disturbance  on  the  free  surface  of  the  magnetic  fluid,  and  determines 
as  well  the  stability  of  the  disturbance  by  considering  the  change  of  type  of 
a  linear  second  order  differential  operator  in  the  equation.  It  has  been 
found  that  in  the  long  wave  limit  the  component  of  the  applied  magnetic  field 
aligned  with  the  direction  of  wave  propagation  has  a  stabilizing  effect  and 
the  transverse  component,  however,  has  a  destabilizing  effect.  Furthermore,  a 
critical  Reynolds  number  has  been  explicitly  obtained.  We  have  also  studied 
the  case  when  the  linear  operator  is  parabolic,  that  is,  the  equation  assumes 
the  form 


‘^t  +  -  “yy  =  °  • 

by  a  coordinate  transformation.  This  equation  subject  to  a  prescribed  initial 
condition  was  solved  by  means  of  the  method  of  characteristics  and  in  general 
the  solution  may  break  down  eventually  [1]. 

We  also  note  that  the  work  on  the  inhomogeneous  Burgers  equation  for  a 
viscous  fluid  flow  down  an  inclined  plane  with  a  prescribed  surface 
disturbance  is  in  good  progress.  We  attempt  to  relax  some  of  the  assumptions 
imposed  on  the  prescribed  disturbance  for  the  existence  of  a  pulse-like 
solution  [2]. 
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